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Abstract 

cn 

^ ' Methods from the geometry of nonholonomic manifolds and Lagran- 

>0 . ge-Finsler spaces are apphed in fractional calculus with Caputo deriva- 

tives and for elaborating models of fractional gravity and fractional 
11^ . Lagrange mechanics. The geometric data for such models are en- 

coded into (fractional) bi~Hamiltonian structures and associated soli- 
tonic hierarchies. The constructions yield horizontal/vertical pairs of 
p^ ' fractional vector sine-Gordon equations and fractional vector niKdV 

equations when the hierarchies for corresponding curve fractional flows 
are described in explicit forms by fractional wave maps and analogs of 
Schrodinger maps. 
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1 Introduction 

The goal of this paper is to show how fractional solitonic hierarchies 
can be canonically generated in various models of fractional gravity and 
geometric Lagrange mechanics. Such constructions are possible in explicit 
form for a class of fractional derivatives resulting in zero for actions on 
constants, for instance, for the Caputo fractional derivative [H [2 [3l [H O [6] . 
This property is crucial for constructing geometric models of theories with 
fractional calculus even, after corresponding nonholonomic deformations, we 
may prefer to work with another type of fractional derivatives. 

This is the second partner work of paper [7] (we also recommend readers 
to consult in advance the papers [H [9l [lOl [11] on details, notation con- 
ventions and bibliography) where we proved an important result that via 
nonholonomic deformations on fractional manifolds and bundle spaces, de- 
termined by a generating fundamental Lagrange/ -Finsler, or an Einstein 
metric, we can construct linear connections with constant coefficient curva- 
ture. For such fractional "covariant" connections, it is possible to provide 
a formal encoding of integer and non-integer gravitational dynamics, Ricci 
flow evolution and constrained Lagrange/Hamilton mechanics into hierar- 
chies of solitonic equations. 

The most important consequence of such geometric studies is that using 
bi-Hamilton models and related solitonic systems we can study analytically 
and numerically, as well to try to construct some analogous mechanical 
systems, with the aim to mimic a nonlinear/fractional nonholonomic dy- 
namics/evolution and even to provide certain schemes of quantization, like 
in the " fractional" Fedosov approach [9l [12] . 

This work is organized in the form: In section [2 we remember the most 
important formulas on Caputo fractional derivatives and nonlinear connec- 
tions. Section [3] is devoted to definition of basic equations for fractional 
curve flows. The Main Theorem on fractional bi-Hamiltonians and solitonic 
hierarchies is formulated and proved in[31 Finally, we derive in general form 
the corresponding nonholonomic fractional solitonic hierarchies in section [5l 

2 Caputo Fractional Derivatives and N— connecti- 
ons 

We summarize some important formulas on fractional calculus for non- 
holonomic manifold elaborated both in global and coordinate free forms, 
as well with important local integro-differential parametrizations, in Refs. 



[TOl \TT\ [HI [7] . Readers are recommended to study in advance those works 
(and references therein) on fractional differential geometry and applications. 
Such a calculus is nonlocal both on space and/or time coordinates when the 
algebra of fractional derivatives does not have the same properties as in the 
integer case. Nevertheless, having well-defined concepts of integral calculus 
for curved nonholonomic manifolds and bundles of integer dimension, we can 
introduce such algebras in local forms and then globalize the constructions 
using corresponding charts on atlases covering corresponding spaces. 

a 

Our geometric arena consists from an abstract fractional manifold V (we 
shall use also the term "fractional space" as an equivalent one enabled with 
certain fundamental geometric structures) with prescribed nonholonomic 
distribution modeling both the fractional calculus and the non-integrable 
dynamics of interactions. We note that in our works a corresponding system 
of notations is elaborated in a form to unify the approaches on fractional 
calculus, nonholonomic bundle spaces, nonlinear connection (N-connection) 
formalism etc. There are considered boldface symbols, over/under and left 
up/low labels etc as we shall explain below. 

Let us consider that f{x) is a derivable function f : [ ix, 2x] -^ M, for 
M 9 Q > 0, and denote the derivative on x as (9a; = d/dx. We note by 
ix and 2X (with left low labels, respectively, 1 and 2) two ends of a real 
line interval. The fractional left, respectively, right Caputo derivatives are 
denoted in the form 



ix 
2X 



X 

For instance, we emphasize that the integral is considered from ix to x 

a 

in the symbol of partial derivative ixd^- We shall always put a over a 
symbol (or up-left/-right to such a symbol) in order to emphasize that the 
constructions are considered for a fractional calculus with a € (0, 1). There 
will be underlined some corresponding symbols if their definition is strictly 
related to the concept of Caputo derivative. Using operators ^ generalized 
on M"', where n = 1,2..., and the same fractional a is associated to any 
such integer dimension, we can construct the fractional absolute differential 

d := [dx^)^ odj when dx^ = (dx^)"' Y(2~a) ' "^^^re we consider ix* = 0. 



We denote a fractional tangent bundle in the form TM for a € (0, 1), 
associated to a manifold M of necessary smooth class and integer dim M = 
n|^ Locally, both the integer and fractional local coordinates are written 
in the form u^ = {x^ ,y"'), where indices to coordinates run values i,j, ... = 
1,2, ...n (for coordinates on base manifold) and a,b, .. = n + l,n + 2,...,n + n 
(for typical fiber coordinates) for integer dimensions but keep in mind that 
the local derivatives are of fractional type ([1]), associated respectively to 

any such integer coordinate. A fractional frame basis ea = e o{u")d_a( on 

TM is connected via a vierlbein transform e Ju'^) with a fractional local 
coordinate basis 



a /a a a '^ \ 

dpi = [dj, = ^^ji dj,,df,i = ^yb' dj,, j , 



for j' = 1,2, ...,n and b' = n-\-l,n + 2, ...,n + n. The fractional co-bases are 
written e = el, {u")du" , where the fractional local coordinate co-basis is 

d'u''' = ('(dx^')°,(d/)°V (3) 

a 

It is possible to define a nonlinear connection (N-connection) N for 

a 

a fractional space V by a nonholonomic distribution (Whitney sum) with 

a a 

conventional h- and v-subspaces, KV and uV, 



a a 



TV = hN®vy. (4) 

Locally, such a fractional N-connection is characterized by its local coeffi- 
cients N={ "TVf}, when N= "iVf (u)(da;^)" oi^. 

a 

On V, it is convenient to work with N-adapted fractional (co) frames, 



°e,=|- -iV;i, °e, = gfe 



(5) 



"e^ = ["e^' = (dx-')", V = ((iy^)"+ "iVfc^(dx'=)"]. (6) 

A fractional metric structure (d-metric) g = { '^gajs} = [ "^dkj-, "ffcfe] on 



'^The symbol T is underlined in order to emphasize that we shall associate the approach 
to a fractional Caputo derivative. 



V can be represented in different equivalent forms, 

g = ^g,_p{u){dvXr®{d^ilr (7) 

= Vk'j' "e''' (g) "e^'' + r/c'fe' "e'"' (g) °e^', 

where matrices r/^/j/ = ^^^[±1, ±1, ..., ±1] and r/a'f,' = dia5[±l, ±1, ..., ±1], 
for the signature of a "prime" spacetime V, are obtained by frame transforms 

Vk'j' = e^k' ^ j' "fi'fcj and ria'h' = e^^, e^y "^gab- 
ex 
We can adapt geometric objects on V with respect to a given N-connection 

a 

structure N, calling them as distinguished objects (d-objects). For instance, 

a a 

a distinguished connection (d-connection) D on V is defined as a linear con- 
nection preserving under parallel transports the Whitney sum dH). There is 
an associated N-adapted differential 1-form 

"r-^ = "r-^^ "e^ (8) 

parametrizing the coefficients (with respect to ([6]) and ([5])) in the form 



"pT I Qri ara a.t~<% curt a 

*- Tii—\ ^jk^ ^bk^ '-'jc^ '-'bcl- 

The absolute fractional differential "d = ixdx-\- ^yfij, acts on fractional 
differential forms in N-adapted form. This is a fractional distinguished 
operator, d-operator, when the value "d := "e" "e^ splits into exterior h- 
and v-derivatives when 

, J, := {dxT ixk = "e^ "e,- and ,ydy := {dyT i Ja = "e^ "e^- 

Using such differentials, we can compute in explicit form the torsion and cur- 
vature (as fractional two d-forms derived for ([8])) of a fractional d-connection 

D = { "r-^^}, 

a^T ^ D "e^ = "d "e^ + "r^^ A "e'' and (9) 

"7^^^ = D °r^^ = "d T^p - "r^^ a "r^^ = "r^^^^ "e^ a "e^ 

Contracting respectively the indices, we can compute the fractional Ricci 
tensor "TZic = { "R^/^ = "R'^^^^} with components 

a Tj _L. OLTjk art _j_ _ a rtk art _^ a jjb art _^ a rtc ('\r\\ 

ixij — n jjf^, n,ia — -Tt j^,^, rtai — -Tt ^jf,, rtab — J^ abc l^"-"/ 



a 

and the scalar curvature of fractional d-connection D, 

with "g"^^ being the inverse coefficients to a d-metric ([7])q 

The Einstein tensor of any metric compatible D, when D,- "g'^P = 0, is 
defined "£'ns = { "Gq,/3}, where 

"Gai3 ■= °R-a^ - 2 "Sa^ "R- (12) 

a 

The regular fractional mechanics defined by a fractional Lagrangian L 

a Q 

can be equivalently encoded into canonical geometric data ( ^N, ig, "D), 
where we put the label L in order to emphasize that such geometric objects 
are induced by a fractional Lagrangian as we provided in [101 El El E] . We 

a 

also note that it is possible to "arrange" on V such nonholonomic distri- 

butions when a d-connection qD = { q^'^^^iri} is described by constant 
matrix coefficients, see details in [HI [14], for integer dimensions, and [7], for 
fractional dimensions. 

3 Basic Equations for Fractional Curve Flows 

In symbolic, abstract index form, the constructions for nonholonomic 
fractional spaces with correspondingly defined distributions are similar to 



■^For applications in modern gravity and geometric mechanics, we can considered more 
special classes of d-connections: 

• There is a unique canonical metric compatible fractional d-connection "D = 
{ "f^,^ = ( "L)k, "Lt„ "q„ "Cfc",)}, when °D ( "g) = 0, satisfying the con- 
ditions "f)^ = and "fl^ = 0, but "f)^, "'f'^^i and "f^^ are not zero. The 
N-adapted coefficients are explicitly determined by the coefficients of 0, 

^jk = 2 9 \ f^k 9jr+ Bj Qkr- 6^ Qjk), 

Lbk = £(,( iVfe) + - ff ( efc gba- Qdc eb N^ - Qdb e^ Nk), 

a ^i ^ a ik a a. a p^a ^ a, ad / a a. , ct a a ct \ 

• The fractional Levi-Civita connection "V = { °T^^o} can be defined in standard 
from but for the fractional Caputo left derivatives acting on the coefficients of a 
fractional metric ((7|. 



those for the Riemannian symmetric-spaces soldered to Klein geometry of 
"integer" dimension. The fractional structure is encoded into the local Ca- 
puto derivatives. 

Following the introduced Cartan-Killing parametrizations, we analyze 
the flow 7(t, 1) of a non-stretching curve in Vn = G/50(n)© SO{m) 

a [-1 

extended to Vn = G/50(n)© 50(m)o This extension is defined via a 
coframe °e G T* Vn O (/ip©fp), which is a N-adapted {SO{n)®SO{m))- 
parallel basis along 7, and its associated canonical d-connection 1-form 
"^r G T* Vn (^ (so(n)©so(?7i)). Such d-objects are parametrized: 



"ex= "e;,x+ "e^x, 
where (for (1, 7^) G M", "o^ G M"~^ and (1, tf) G I 



and 



"e^x = 7hxJ/i "e 



'et,x = 71'xjt' °e 







\ G 
/lO 



»m— 1^ 



(1,0: 



^We use an isomorphism between the real space so{n) and the Lie algebra of n x n skew- 
symmetric matrices. This allows us to establish an isomorphism between hp ~ R" and 
the tangent spaces T^M = so(n + 1)/ 5o(n) of the Riemannian manifold M = SO(n + 1)/ 
SO{n) as described by the following canonical decomposition 



hg — so{n + 1) Z) hp e 










for /lO £hi} — so(n) 



with /ip = {p' } SR" being the h-component of the d-vector p = (p' ,p" ) and hp^ mean 
the transposition of the row hp. In our approach, TxM — >■ T_^M, with Caputo fractional 
derivatives. The Cartan-Killing inner product on /ig is 



hp-hp = 




-/tp^ 



hp 
hO 




-/ip^ 



-tr 




-hp^ 



hp 
hO 



hp 
hO 


-hp'^ 



hp 
hO 



where tr denotes the trace of product of matrices. This product identifies canonically hp 
~ R" with its dual hp* ~ R". In a similar form, we can consider 



VQ = so{m + 1) 3 i;p G 




-up 



vp 
vO 



for vO Evt) = so(7Ti) 



with vp — {p" } £R™ being the v-component of the d-vector p — (p* ,p'^ ) and define the 
Cartan-Killing inner product wp-iip =|ir{...}. In general, we can consider the Cartan- 
Killing N-adapted inner product p ■ p =hp-hp+vp-vp. 



We also introduce 

for 

r/ix=7fexjL 



°r = [r/ix,r^x] , 





(0,^) 



eso(n + 1) 



where L 



if 



-(0,0^)^ L 
e so(n), if e IR"'"^ hO eso(n - 1), and 



v^K.'- 



=7.xJC 





-(o,tr 



c 



where C 



t7T ^0 







G so(m), V G 



pm— 1 



G 50(171+ 1), 



t;0 Gso(m — 1). 



There are decompositions of horizontal SO{n+ 1)/ SO{n) matrices, 



hp 3 



+ 



hp 

-hp^ m 


-(0,/iF^_ 



/ipii, 



/ip||, 



\T 



(0,/i^_ 

m 



into tangential and normal parts relative to "e/^x via corresponding decom- 
positions of h-vectors hp = (/ip||,/ip_L) gM" relative to I 1, I , when hp\\ 

is identified with hpc and /i p _l is identified with /ip_|_ = /ip^± . In a similar 
form, it is possible to decompose vertical SO{m + 1)/ SO{m) matrices. 



vp B 





vp 

-vp'^ vO 


= 




. -(^Pll'"^) 


(vp\\, 


f 





)" 


t;0 


? 





into tangential and normal parts relative to "e^x via corresponding de- 
compositions of h-vectors vp = (t;p||,t;p^) gM™ relative to I 1, I , when 

t;p|| is identified with vpc and vp± is identified with vp± = vpQ±. Locally, 
we consider, for instance, instead of M" the fractional space "M" local (co) 
vectors defined by Caputo fractional derivatives and their duals. 



The canonical d-connection "D = { '^F'^o } [7] induces matrices decom- 
posed with respect to the fractional flow direction. In the h-direction, we 
parametrize 



*e/jY = 7tJ/i "e 



{he\\,h'^±) 



-{he\\,h'^_i) hO 

when "e/jY e hp, {he\\,h'tj_) £ "M" and h'^± G "M"-i, and 



r/iY=7/iYjL 



(0,"^) 

-(0,1^)'^ hvJr 



where hzu-j-- 



^ 



G 50{n+ 1), 
Gso(n), ^G °M"-i, /leGso(n-l). 



(13) 



In the v-direction, we have 
e„Y =7rJ^^e = 



tr \^ 



(feiijV e _|_) 



(veil, 11 e _l) 



t;0 



when e^Y € vp, (veu,v^±) G "M™ and u^^ G "R"'"^ and 



rt,Y=7»;YjC 



(0, 



G so(m + 1) 



where vzur- 



fe 



Gso(m), feG °M"^-\ weGso(m-l). 
The components heii and he ± correspond to the decomposition 

into tangential and normal parts relative to "e/jx- In a similar form, one 
considers t;e|| and v e ± corresponding to the decomposition 



*e^Y = fg(7r;7i) °e^x + (7r)-Lj^e^- 



Working with such matrix parametrizations, we define 



[ "e;,x, "e;,Y] = - 


- 


" ■ 

/lei 


E so(n+ 1), 






(] 


for he± = 


-1 

h'^i_ ' 


Gso(n); 








[ThY, "e;,Y] = - 


- 


(0,^)" 
_-(0,^f 


e/ip±; 








[r;.x, "e^Y] = - 


- 


(- 

-{-it ■h'^^,he\\l^f 


-It 


hO 


/ie||l^) " 


£hp- 








and 










[ "e^x, "e^v] = - 


- 


" " 

t;e_L 


€ 50{m + 1), 









for ve± = 


v'^± ' 


e so(m); 








[ri,Y, "g^jy] = - 


- 


(0,fe)" 
. -(0,fef 


evp±; 








[r.x, "e,Y] = - 


- 


_-(-V. 




V e ^,v 


e||t 


,'- 


-t 


vO 


uey V) 



e vp. 



We use formulas (\T4^ and p5]) in order to write the structure equations in 
terms of N-adapted curve fractional flow operators soldered to the geometry 
Klein N-anholonomic spaces, the formulas are "fractional" extensions of 
the respective ones in Refs. [131 El- This way, it is possible to construct 
respectively the G-invariant N-adapted torsion and curvature generated by 
the fractional canonical d-connection, 



'T(7.,7i) 



^Dx7r- "Dy7i J °e 



OTI a 



Dx "eY - "Dy "ex + 



rx, "eY 



(16) 
^fv, "ex 



and 



^R(7.,7i)"e 



Dx, "Dy 



(17) 



a-r\ Oi 



Dx "Fy - "Dy "rx + 



o-i-i a 



rx, "Fy 



10 



where "ex = 7iJ "e, "ev = 7rJ "e, "Fx = 7iJ "F and "Fy ^7rJ °r. 

Applying a d-connection "D (in particular, we can take °D) instead 
of the Levi-Civita one °V, we get 

= ("D;,x7r- °D;,Y7i)J/i°e (18) 

= "D/ix "e/iY — "DftY "e/jx + [ "L/^x, "g/iy] - [ '^^hY, °e/,x] ; 
= ("D„x7r- "D„Y7i)J^'"e 

= "Dt)X "e^jY — "^D^Y ""e^x + [ °C„x, ""g^jy] — [ °C„y, "e„x] , 

/i "R(7r,7i)/i "e = [ °D;,x, °D;,y] h "e 
= "D/iX "LftY — "D/iY "L/iX + [ "L/iX, "L/iy] 
^ "R(7r,7i)^^ "e = [ "D,x, "D,y] ^ "e 

= -U^x ^vY — L>vY *^v:X. + ^vJi., ^vY\ ■ 



Following N-adapted curve flow parametrizations (J14p and (jlSp . the 
equations ([TH]) are written 



0= °D/,x/ie|| +lf -/i"^^, 0= °D^,xwe|| + V-v1e"_L; (19) 

= ^ - /leylt + "D/jx/i'^±, = fe - we||V + "D„x^^^±; 
"D,,x^- "DfeY^ + ^J/i0 = /i'^±, 

For such fractional spaces, the tensor and interior products, for instance, for 
the h-components, are defined in the form: (8) denotes the outer product of 
pairs of vectors (1 x n row matrices), producing n x n matrices A B = 
A B, and J denotes multiplication of n x n matrices on vectors (1 x n 
row matrices) ; one holds the properties A\lB(S>C] = iA-B]C which 
is the transpose of the standard matrix product on column vectors, and 
(b^C) A = (C- A] B. As basic vectors, we use the Caputo fractional 
derivatives. Similar formulas hold for the v-components but, for instance, we 
have to change, correspondingly, n — > ?n, and A ^ A; for such constructions 
the fractional differentials have to be used. 

Lemma 3.1 On nonholonomic fractional manifolds with constant curva- 
ture matrix coefficients for a d-connection, there are N-adapted fractional 
Hamiltonian symplectic operators, 

h^J= °Dfcx+ "D-^ (It-) It and v^J= "D,x+ "D;^ (V-) V, (20) 

11 



and cosymplectic operators 

h'"H= °Dftx + l^J "D^^ (l^A) and v '^U = °D,x + Vj "D;-^ ( Va) , 

(21) 
where, for instance, A A B = A (g) B — B® A. 

Proof. We sketch some key steps of the proof. The variables en and 0, 
written in h- and v-components, can be expressed correspondingly in terms 
of variables v,w,he± and v ,w,v e ± (see equations (fT9|) ). 

tell = - "D-^(lt • h-^±), t;e|| = - "D;;^(V • 7;iFx), 

and h& = "D"^ (it O ^ - ^ O it) , u© = "D~^ (V O fe - fe O V) . Sub- 
stituting these values, respectively, in equations in (|19p . we express 

fe = - "D,x^^± - "D;x(^ • ^^±)^, 

contained in the h- and v-flow equations respectively on v and v , consid- 
ered as scalar components when "^DiiyV = Vr and "D/^y v = v r, 

it^ = "D,,x^-ltj °D--^(lt0^-^«)lt)- "^/i-^x, (22) 
V^ = "D^xfe- Vj "D-^(V«)fe-fe«)tr)- "^'(;iS"^, 

where " R and " S are the scalar curvatures of chosen d-connection. For 
symmetric Riemannian spaces like SO{n + 1)/S0{n) ~ S", the value R is 
just the scalar curvature x = 1- Oii N-anholonomic fractional manifolds, it 
is possible to define such d~connections that " R and " S are certain zero 
or nonzero constants. D 

The properties of operators (f20|) and ([2T]) are defined by 



Theorem 3.1 The fractional d-operators °'J = {h °'J,v °'J) and 
"'Ti = {h '^'H,v "7i) are respectively {0{n — l),0{m — 1)) -invariant Haniil- 
tonian symplectic and cosymplectic d-operators with respect to the frac- 
tional Hamiltonian d-variables (17, v ) . This class of d-operators defines the 
Hamiltonian form for the curve fractional flow equations on N-anholonomic 
fractional manifolds with constant d~connection curvature: the fractional h- 
flows are given by 

it, = /i"?{(^)- "^/i^x = /i"9^(/i'^x) - "^^ /i"^^, 
^ = /i°j(/i-^^); (23) 

12 



the fractional v-flows are given by 

fe = v'^Jivt^), (24) 

where the so-called (fractional) heriditary recursion d-operator has the re- 
spective h- and v-components 

h''^ = Wnoh'^J and v''^ = v''nov "J. (25) 

Proof. Such a proof follows from the Lemma and (|22|) . D 

4 Fractional Bi— Hamiltonians and Solitonic Hier- 
archies 

The fractional recursion h-operator from (|25|) . 

h'^m = "D,,x ( "D;,x + "D^i (t-) t) + If J "D^i (If A "D;,x) 

= "D|x + I "D,xP + "D,i (If.) If, - If J "D-i(lfiA), (26) 

induces a horizontal hierarchy of commuting Hamiltonian vector fields he^ 

starting from he j_ = v\. Such vector fields are given by the infinitesimal 
generator of 1-translations in terms of arclength 1 along the curve. A vertical 
hierarchy of commuting vector fields v e j_ starting from v e j_ = u i is 
generated by the recursion v-operator 

V ^m = "D,x ( "D,x + "D;^ (V-) V) + Vj "D;^ (V a "D,x) 

= "D^x + I "D.xP + "D,-^ (V.) V, - Vj "D;^(V,A). (27) 

We can associate fractional hierarchies generated by adjoint operators "9^* = 
(h °9^*, V °9^*), of involuntive fractional Hamiltonian h-covector fields 
'^(^) = S {h "'H^ >^ jSv in terms of fractional Hamiltonians h °^H = 
h °i/('=)(lf, if 1,^21,-) Starting from ^W = ^ ,h "F^") = i|lf|2 and of 
involutive fractional Hamiltonian v-covector fields w^ ' = 5 {y °H^ >^ / 6 v 
in terms of Hamiltonians v "'H = v °'H^'''{v , v \, V2\,---) starting from 
^(0) = y ^y "^(0) = ■^\ V p. The relations between fractional hierarchies is 
given by formulas 

13 



where k = 0,1,2,.... All hierarchies (horizontal, vertical and their adjoint 
ones) have a typical niKdV scaling symmetry, for instance, 1 — )• Al and v -^ 
X~^v under which the values he j_ and h °^H^ ' have scaling weight 2 + 2k, 
while w^ ' has scaling weight 1 + 2k. 

Corollary 4.1 There are N-adapted fractional hierarchies of distinguished 
horizontal and vertical commuting bi~Hamiltonian fractional flows, corre- 
spondingly, on V and v, associated to the recursion d-operator [25\) given 
by 0{n — 1) (B 0{m — 1) -invariant d-vector evolution equations, 

It. = /i-^f +1) - -t h^f = h^n{5 {h °i7(^-.^)) 15-^ 
= {h ""jy^ U [h "i7(^+i.^)') /5^ 

with horizontal fractional Hamiltonians 

j^ ajjik+l,t) ^ ^ ajj{k+l,t) _ ^th "//(^•'^) 

and 

tr = t;^f +'^ - "^ t;^f =v'''h(^6 (v "//C^.^)) /6 

= {V "^jy^ (b (v ^H^k+l,'s)\ /^^\ 

with vertical fractional Hamiltonians 

for k = 0,1,2, The fractional d-operators "Ti and "J are N-adapted 

and mutually compatible from which one can be constructed an alternative 
(explicit) fractional Hamilton d-operator ""H = ""H o "^ o ""H = "9^o ""H. 

Proof. It follows from above presented considerations. D 

4.1 Formulation of the Main Theorem 

Our goal is to prove that the geometric data for any fractional metric 
(in a model of fractional gravity or geometric mechanics) naturally define a 
N-adapted fractional bi-Hamiltonian flow hierarchy inducing anholonomic 
fractional solitonic configurations. 
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Theorem 4.1 For any N-anholonomic fractional manifold with prescribed 
fractional d-metric structure, there is a hierarchy of bi-Hamiltonian N- 
adapted fractional flows of curves 7(r, 1) = /i7(r, 1) + V7(t, 1) described by 
geometric nonholonomic fractional map equations. The fractional flows 
are defined as convective (traveling wave) maps 

7r = 71, distinguished {h'y)^ = {h'y)f^-^ and {vj)^ = (v7)t,x • (^8) 

There are fractional +1 flows defined as non-stretching mKdV maps 

-ihj)^ = "DL(/i7),x + ^rD,x(/i7).xlL (^7)^x, (29) 
-{v-fl = "D^xM.x + ^rO.xM.xl'g (^7).x, 

and fractional +2,... flows as higher order analogs. Finally, the fractional 
-1 flows are defined by the kernels of recursion fractional operators i26\) and 
Iji27\ ) inducing non-stretching fractional maps 

"D;,Y (/i7)/,x = «^^ "D„Y (t'7).x = 0. (30) 

Proof. It is given below in section I32J D 

4.2 Proof of the Main Theorem 

We generalize for fractional spaces a similar proof from Refs. |131 [T^ 
sketching the key steps for horizontal flows. The vertical constructions are 
similar but with respective changing of h- variables / objects into v- vari- 
ables/ objects. By corresponding nonholonomic constraints we can empha- 
size certain h- and v-evolutions which are inter-related. 

We get a fractional vector mKdV equation up to a convective term (which 
can be absorbed by redefinition of coordinates) defining the fractional +1 
fiow for /i e _|_ = vi, 

when the fractional +{k + 1) flow gives a vector mKdV equation of higher 
order 3 + 2A; on u and there is a h-flow v-r = v\ arising from he ± = 
and he \\ = 1 belonging outside the hierarchy generated by h "IH. Such 
fractional flows correspond to N-adapted horizontal motions of the curve 
7(t, 1) = /i7(r, 1) + vjir, 1), given by 

(/i7). = /((/i7),x, "D,x(/i7),x, "Dk(/i7),x,-) 
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subject to the non-stretching condition | (/i7)/ix Ug ~ 1' when the equation 
of fractional motion is to be derived from the identifications 

(/i7), ^^ "e;,Y, °D;,x (/i7)/,x ^^ "^/.x "e^x = [ ''Ux, "e/,x] 

and so on, which maps the constructions from the tangent fractional space 
of the curve to the space hp. For such identifications, we have 



[ °L/iX) "e/ix] 
[ "L/iX, [ "L/iX, "e/ix]] 







(0,1^) 



(0,lt) 



hO 



ehp, 



-(|lt|2,'0^) /lO 

and so on, see similar calculus in (J14p . Stating for the fractional +1 h-flow 



h~^± = 1^1 and h-^\\ = - "D^-^ (1^ • l^i) = -^I^P, 



we compute 



e/iY 



(/ie||, /i e ^)' 



-lii?p 



[he\\,h'S 



0, 



T 



+ 



(0,Ax) 



1 



= °D?iX [ "L/iX, "e/ix] + - [ °L/iX, [ °L/jx, "e/ix]] 

= - "P,x [ "L;,x, "e^x] - ^I^P "e,x. 

The above identifications are related to the first and second terms, when 

|1^P = < [ "L;,x, °e;,x] , 

[ "L^x, "e^x] >hp^^ h "g ( "D;,x (/i7)/.x ' "D/^x (/i7)/.x) 
= rD,x(/i7),xl^^g> 

and allow us to identify "I'/ix [ "L/ix, "e/ix] to "D^x (^7)ftx • ^^ ^ result, 
we have 

-"e,Y^^ "Dk(/i7)/.x + ^rD/.x(/i7)/.xl^<.g (/^7),x 
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which is just the fractional equation (j29p in the Theorem 14. 1 1 defining a non- 
stretching mKdV map h-equation induced by the h-part of the fractional 
canonical d-connection. 

To derive the higher order terms of hierarchies, we use the adjoint 
representation ad{-) acting in the Lie algebra hg = hp ©so(n), with 



where v 
we get 










0, ot 



Gso(n + 1) 







e so n 



Applying again a(i(["L,,x, "e/,x]) , 



ad(["LfcX, "ehx])' "e;,x 



1^12 a 



1, 



T 



1,-cf 



e/ix, 



when the fractional equation (I29p can be represented in alternative form 



(^7). 



"D|x (/i7) 



feX 



^"i^-iad("D,x(M,x)' (M 



hX 



Finally, we consider a fractional -1 flow contained in the h~hierarchy de- 
rived from the property that h e ±is annihilated by the fractional h-operator 
h "'J and mapped into h "9^(/i e _|_) = 0. This mean that h °'J{h e _|_) = 
ro = 0. Such properties together with (fT3|) and fractional equations (p2]) 
imply "^Lt- = and hence 

h "P. "e;,x = [ "L,, °e/,x] = for /i "P, = /i "D, + [ °L„-]. 

We obtain the equation of fractional motion for the h-component of curve, 
/i7(r, 1), following the correspondences "D/jy < — > h "^Dr and /171 < — > 
"g/iXi "D/jy (^7(''") 1)) = 0, which is just the first fractional equation in 
(I30D. 



5 Fractional Nonholonomic mKdV and SG Hier- 
archies 

In this section, we consider explicit constructions when fractional soli- 
tonic hierarchies are derived following the conditions of Theorem 14.11 
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The fractional h-flow and v-flow equations resulting from ()30p are 

^r = - "^/i'^± and V^ = - "^Tji^i, (31) 

when, respectively, 

= ^ = - "D,,x/i'^± + /ie||lt, "D/jx/ie|| = /i"^± • it 

and 

= fe" = — "D^xf^ e ^ + t;e|| V, "Dt,xwe|| = u e _|_ • t7. 

The fractional d-flow equations possess horizontal and vertical conservation 
laws 

"D,x((/ie||)2 + [/i^^|2)=o, 

for (/ie||)2 + |/i'5"±P =< her, her >hp= I (^7)r ll ^^g' and 

°D„Y((we||)2 + |^1S-^|2)=o, 
for (1^611)^ + iv e _|_p =< ver,ver >vp= \ (vj)^ \l a„. This corresponds to 

A |2 _g 

^t I (h-^) I? . 

g 



'D;,x|(/i7)rlh "K = and "D,x|M,l' ^.g = 0. 



We can also rescale conformally the variable r in order to get | {hj)^ || ^ 
1 and (it could be for other rescaling) | (^7)^ \l a™ = 1, i.e. to have 



|2 
It) 

\2 , \u^ |2 _ 1 „„, / .„ \2 , i„.tr |2 



{he\\y + \h'g ±\^ = 1 and (vey)^ + jt;"© ^|^ = 1. 

Then, we express heii and /i e _|_ in terms of v and its derivatives and, sim- 
ilarly, we express ve\\ and v e _l in terms of V and its derivatives, which 
follows from (|3ip . The N-adapted fractional wave map equations describ- 
ing the -1 flows reduce to a system of two independent nonlocal fractional 
evolution equations for the h- and v-components, 



"vr = - "D,;;i a/ «^2_|V J2t7 



We can rescale the equations on r to the case when the terms ° it^, " 5 ^ = 
1, and the fractional evolution equations transform into a system of hyper- 
bolic d-vector equations, 



^D;,x(ltr) = "Vl - l^rP ^ and "D,x(V,) = -^1 " I^^P V, 

(32) 
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where "D^x = °^dh\ and "D^,x = "'dyi are usual partial derivatives on 
direction 1 =hl+vl with Vr and V,- considered as scalar functions for the 
covariant derivatives "D/jx and "D„x defined by the fractional canon- 
ical d-connection. It also follows that he ± and v e ± obey corresponding 
fractional vector sine-Gordon (SG) equations 

U{i-\h^^\2yi »dHi{h^±)] = -h-^i. (33) 

l_|„i^^|2)-i ''d^iiv^^)] = -v^^. (34) 



Conclusion 5.1 The recursion fractional d-operator "IH = (/i "D\,h °1H) 
i25\) . see i26\) and \21l^ , generates two hierarchies of fractional vector rriKdV 
symmetries: the first one is horizontal, 

lf(0) = ltM,^^r^=h ''m{lthl) = 1^3hl+l\^\'^M, (35) 

+ 1 ((I^P)m hi + |lt.iP + ^11^1^) Ai - ^1 Alp If, 

with all such terms commuting with the fractional -1 flow 

{ifry^ = h'^^ (36) 



associated to the fractional vector SG equation ^3^): the second one is ver- 
tical, 



V(0) 



3, 



V,i, V« = V "^(V,i) = V3.1 + -|Vp V,i, (37) 



2' 



5 



+ 1 ( (I^Pk ^'l + I Alp + tI^P ) Al - ^1 Alp V, 



2 V ' ' ' ' 4' 



') ^-41 



wit/i a// sitc/i terms commuting with the fractional -1 flow 

(A)-' = v"^^ (38) 

associated to the fractional vector SG equation \°J4-^ . 
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Proof. It follows from the above, in this section, and Corollary 14. li D 
Finally, using the above Conclusion, we derive that the adjoint frac- 
tional d-operator "IH* = °^ J o "7^ generates a horizontal hierarchy of 
fractional Hamiltonians, 

h "ifw = i|^|^ h "//« = -^|lt,lp + i|l^|^ (39) 

and vertical hierarchy of fractional Hamiltonians 

^; -//(O) = 1|V|2, z; "i/W = -i|V.,p + i|V|^ (40) 

2 2 o 



z 4 z io 
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all of which are conserved densities for respective horizontal and vertical 
fractional -1 flows and determining higher conservation laws for the corre- 
sponding hyperbolic fractional equations (f33|) and (p4]l . 

Finally we note that two concrete examples of fractional solitonic so- 
lutions in gravity are studied in section 5.3 of Ref. [TT] . The "integer" 
solitonic hierarchies from |13l [H] can be similarly generated in " fractional" 
forms by using corresponding fractional Caputo derivative operators. The 
length of this paper does not allow us consider such applications. 

Acknowledgement: This work contains the results of a talk at the 
3d Conference on "Nonlinear Science and Complexity", 28-31 July, 2010, 
Qhankaya University, Ankara, Turkey. 
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